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Abstract 

We review results obtained with a recently proposed variational cluster approach (VCA) for the competition 
between d-wave superconductivity (dSC) and antiferromagnetism (AF) in the high-Tc cuprates. Comparing the 
single-particle spectra of a two-dimensional Hubbard model with quantum Monte-Carlo (QMC) and experimental 
data, we verify that the VCA correctly treats the low-energy excitations. The cluster calculations reproduce the 
overall ground-state phase diagram of the high-temperature superconductors both for electron- and hole-doping. In 
particular, they include salient features such as the enhanced robustness of the AF state in case of electron doping. 
For electron- but also for hole-doping, we clearly identify a tendency to phase separation into a mixed AF-dSC 
phase at low and a pure dSC-phase at high doping. 
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The central issue in the field of high-temperature su- 
perconductivity (HTSC) is the connection of the mi- 
croscopic interactions at the level of electrons and ions, 
which are at high energy and temperature T, with 
the "emerging phenomena" at T = 0, i.e. compet- 
ing and nearly degenerate orders - antiferromagnetism 
(AF), d-wave superconductivity (dSC), heterogeneous 
phases, etc. We will not go into a lengthy discussion of 
what interactions should be retained at the electron- 
ion level. But, when choosing the two-dimensional (2D) 
one-band Hubbard model, i.e. 

where ti t j denote hopping matrix elements, rij,T the 
density at site i with spin ""f" and U the local Coulomb 
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repulsion, one has introduced gross simplifications, 
leaving out other orbital (e.g. p) degrees of freedom, 
long-range Coulomb interaction, electron-phonon cou- 
pling, etc [1]. Nevertheless, this model choice appears 
to be legitimate, last not least in view of the amazing 
agreement achieved between numerical simulations 
and experimental results for the normal-state proper- 
ties of the cuprates [2,3]. 

At low temperatures different orders appear, which 
are not separated by distinct energy scales but com- 
pete with each other. What is required is a kind of 
"magnifying lens" which allows to resolve these com- 
peting orders. Ideally, one should employ a systematic 
renormalization-group approach to integrate out the 
irrelevant degrees of freedom and, thereby, to correctly 
bridge high to low energies and eventually to go to 
T = 0. For the strong-correlation case, realized in the 
HTSC, how to do this is, however, by no means obvi- 
ous. In this context, cluster techniques provide an al- 
ternative way to systematically approach the infinite- 
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size (and, thereby, low-energy) limit. 

Here, we review progress obtained with the varia- 
tional cluster approach (VCA), which was proposed 
and used by Potthoff et. al. [4,5]. This approach pro- 
vides a rather general and controlled way to go to the 
infinite-sized lattice fermion system at low tempera- 
tures and at T = 0, in particular. The ground-state 
phase diagram of the 2D one-band Hubbard model was 
calculated within VCA by Senechal et. al. [6] and, in- 
dependently, by two of us [7]. There are certain tech- 
nical differences, which we discuss below, but the "up- 
shot" of the two works is as follows: For the cluster 
sizes used in the VCA, the T — phase diagram of the 
Hubbard model (1), with hopping terms up to third- 
nearest neighbors, correctly reproduces salient features 
of the HTSC, such as the AF and dSC ground states 
in doping ranges, which are qualitatively in agreement 
with electron- and hole-doped cuprates. 

The VCA is based on the self-energy-functional 
approach (SFA) [8]. The SFA provides a variational 
scheme to use dynamical information from an exactly 
solvable "reference system" (for example an isolated 
cluster) to approximate the physics of a system in the 
thermodynamic limit. For a system with Hamiltonian 
H — Ho(t) + Hi (17) and one-particle and interaction 
parameters t and U , the grand potential is written as 
a functional of the self-energy E as 



M 



fit,C7[£] = Fu[S] + Trln (Gq] - £) 



(2) 



with the stationary property <5fit,t/[£phys] = for the 
physical self-energy. Here, Go,t = (w + — t) _1 is the 
free Green's function at frequency ui, and \i is the chem- 
ical potential. Fxj [£] is the Legendre transform of the 
Luttinger-Ward functional and determines the fully in- 
teracting Green's function via G = — SFu [£]/<$£. It is 
important to note that Fu [£] is a universal functional: 
The functional dependence is only determined by the 
interaction parameters U (for example, the Hubbard 
interaction in Eq. (1)). Therefore, the functional Fu [£] 
is the same as the functional for a problem which is 
"simpler" and solvable, i.e. for a Hamiltonian H' — 
Ho(t') + Hi (17) with the same interaction part but a 
one-particle part that makes it exactly solvable. The 
stationary solutions are obtained (and this is the ap- 
proximation) within the subspace of self-energies E = 
E(t') of that simpler solvable problem that is spanned 
by varying t' . If one takes a single site and connects 
it to continuous (non-interacting) bath degrees of free- 
dom (another H' choice), one recovers the dynamical 
mean-field theory (DMFT) [9] or, for a cluster of sites 
connected to a bath, a cluster variant of DMFT [10]. 

In the VCA, considered in the following discourse, 
H' is build up of disconnected clusters, which have 
their inter-cluster hopping terms removed. In our T = 
approach, the isolated cluster is solved by exact di- 
agonalization. Its Hamiltonian H' includes additional 
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Fig. 1. Upper part: Density plot of the spectral function for 
the 2D Hubbard model at half-filling, T — and U = St 
(t: nearest-neighbor hopping) as obtained by the VCA [5]. 
The lattice is covered by VTO X vTO clusters. Bottom: QMC 
(maximum entropy) result, taken from Ref. [3], for the same 
parameters but for a finite low temperature T — O.lt and an 
isolated 8x8 cluster. Dark (light) areas correspond to large 
(small) spectral weight. 

symmetry-breaking "Weiss" fields [5] to allow for long- 
range order. The VCA solution is finally obtained as a 
stationary point determined by dQt,u[^{t')]/dt' — 0. 
Is this a controlled route to a (T = 0) infinite-size ap- 
proach? To answer this, consider a few "tests" : 

(i) The VCA correctly reproduces long-range AF or- 
der in 2D and the absence of this order in ID [5]. This 
non-trivial test implies that the VCA goes well beyond 
ordinary mean-field theory. 

(ii) An advantage, compared to variational schemes 
based on wave functions [11] , is that the VCA quite nat- 
urally gives the one-particle Green's function G. Fig. 1 
compares the spectral function A(k, uS) oc ImG(fe,w) 
of the VCA for the 2D Hubbard model at U = 8f, half- 
filling and T — with corresponding low-temperature 
QMC data [3] for an isolated 8x8 cluster. One can 
clearly see that the VCA, with the lattice covered by 
\/T0 x \/T0 clusters, correctly reproduces coherent and 
incoherent "bands" (known from ARPES data [12]). In 
particular, the non-trivial proliferation of AF spin cor- 
relations from cluster to cluster, which builds up the 
coherent quasi-particle "band" is obviously correctly 
embedded in the VCA [5]. Similar calculations have 
been performed for the spectral function A(k, ui) of the 
hole- and electron-doped Hubbard model [6,7]. The 
characteristically different doping dependencies give 
rise to different Fermi-surface evolutions upon doping. 
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Fig. 2. Antiferromagnctic and superconducting order param- 
eters, m and A, and chemical potential as functions of hole 
doping x. A and are plotted for the AF+SC (green, Aaf+SC; 
/^af+sc) as well as for the pure SC homogeneous solutions 
(blue, Asc, Msc)- Note that A is scaled by a factor 10 for 
convenience. For x < x\ the system exhibits a coexistence of 
AF and dSC order. Phase separation (PS) occurs between the 
doping levels x\ and X2- For x > X2 pure dSC is realized. In 
the phase separation region x± < x < X2, the homogeneous 
solutions become unstable, and the system prefers to separate 
into a mixture of two densities corresponding to Xi and x^. 
The chemical potential [i c is determined by the Maxwell con- 
struction shown in the figure. At the slope of the AF+SC 
solution changes sign. 



8t and next-nearest-neighbor hopping t n . n .n = — 0.3i, 
obtained with our calculation, is plotted in Fig. 2 for 
the hole-doped and in Fig. 3 for the electron-doped 
case. In the upper part of each figure, we display the 
AF (m) and dSC (A) order parameters as a function of 
doping x. In the lower part of the figures, the chemical 
potential fi is plotted as a function of x. 

Let us discuss hole doping first (see Fig. 2). For dop- 
ings x below a critical value x\ we find a homogeneous 
symmetry-broken state in which both, the AF as well 
as the dSC order parameter m and A are non-zero. 
This corresponds to a phase AF+SC where AF and 
dSC order microscopically and coherently coexist. A 
homogeneous state with pure dSC (m = and A > 0) 
is obtained for dopings x > X2- 

Fig. 2 also shows A and /i for the homogeneous 
AF+SC and SC phases in the range x\ < x < X2 - Here, 
however, these phases are thermodynamically unsta- 
ble. For dopings x with x\ < x < X2, macroscopic 
phase separation between the two phases occurs. In 
practice, doping-dependencies are calculated by vary- 
ing the chemical potential jj,. Following up the grand 
potentials for the two homogeneous phases as func- 
tions of i.e. J7af+sc(m) and Qsc{n), it is found (see 
Ref. [7] for details) that there is a crossing at a critical 
chemical potential fi — \i c (at this point the AF order 
parameter m is still nonzero). Thus, the transition is 
first order as a function of \x. At the transition point 
fi c , the dopings corresponding to the AF+SC and to 
the SC phase, xaf+sc and a?sc, are different. Conse- 
quently, there is a jump Aa; = zaf+sc — £sc at 



Furthermore, the single-particle excitations provide in- 
sight into the characteristic differences in the ground- 
state phase diagram for hole- and electron-doping [7] . 

(iii) To test the stability of the homogeneous phases 
with respect to phase separation (PS), we consider a 
reference system H' of isolated 2x2 clusters where, in 
addition to the two symmetry-breaking terms (Weiss 
fields) H' AF and Hg C , a term H[ ocal is optimized within 
the variational procedure [7]. H' loca[ = £*}2 i(j ni a de- 
scribes a shift e of the chemical potential in the clus- 
ter with respect to the physical chemical potential fi. 
The use of the additional variational parameter e is re- 
quired in order to have a consistent treatment of the 
particle density. The optimization of e has to be done 
simultaneously with the optimization of the parame- 
ters /iaf and hsc, namely the staggered magnetic field 
in the term H' AF and the nearest-neighbor d-wave pair- 
ing field hsc in the term H' sc . Notice that the Weiss 
fields /iaf and hsc are different from the correspond- 
ing order parameters m and A plotted in Figs. 2 and 
3. Quite generally, however, a nonvanishing stationary 
value for the Weiss fields produces a nonvanishing order 
parameter, although the latter can be much smaller. 

The phase diagram for the Hubbard model with U = 
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Fig. 3. The same as Fig. 2 but for electron doping. Note the 
enhanced robustness of the AF state and the strongly reduced 
scale A/x = (fi* — /j, c ) as compared to hole doping. 
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indicating phase separation between a weakly doped 
AF+SC and a higher doped SC phase. 

Due to the inclusion of H[ oca y, \jl c can equivalently 
be obtained by a Maxwell construction. This is shown 
in Fig. 2 where, in the lower panel, \i is plotted as a 
function of x. Here, phase separation is signaled by 
the fact that the j-i(x) is not a monotonous function. 
The Maxwell construction shown in the figure then 
identifies the two dopings x\ and X2 into which the 
system tends to phase separate, as well as the chemical 
potential /j, c in the phase-separated state. In Fig. 2, (i* 
is the point where the slope of n(x) changes sign. For 
H < the AF+SC solution ceases to exist 

Let us now discuss the electron-doped case. While 
the phase diagrams in Figs. 2 and 3 are qualitatively 
similar, the phase in which long-range AF order is re- 
alized is spreading to significantly larger doping values 
in the electron-doped case, in overall agreement with 
the experimental situation. Another important differ- 
ence concerns the energy scale for phase separation, i.e. 
A/i = (/.t* — fj, a ). As one can see from the comparison 
between Figs. 2 and 3, An is an order of magnitude 
larger in the hole-doped case. In Ref. [7] it is argued 
that this can explain the different pseudogap and SC 
transition scales in hole- and electron-doped materials. 
This may give support to theories [13] which are based 
on the idea that fluctuations of the competing phases, 
or of the related order parameters, are responsible for 
the pseudogap phenomenon. 

In order to resolve the relevant small energy scale, 
it is necessary to evaluate Q as well as its stationary 
points with high accuracy. Furthermore, the inclusion 
of the chemical potential shift term H' local consider- 
ably complicates the variational optimization. For the 
rather small clusters of size 2x2 used here, the ref- 
erence system can be treated by full diagonalization 
and the frequency integrals, which are implicit in Eq. 
(2) [8], can be carried out by means of a sum over the 
negative poles of the Green's functions. Senechal et. 
al. [6] have considered clusters up to 10 sites and report 
similar results for x « but, without the inclusion of 
H'i oca i, one cannot reliably test the stability of the ho- 
mogeneous phases against phase separation. Note that 
the transition from the AF+SC to the SC phase may 
appear as continuous as a function of x if phase sepa- 
ration is not taken into account. 

In conclusion, there has been substantial recent 
progress in relating the "high-energy" physics of the 
Hubbard model and its variants to the low-energy 
physics of the competing phases AF, dSC, charge 
inhomogenities, etc. This progress is due to the devel- 
opment of quantum-cluster theories, such as the VCA 
discussed here but also due to cluster extensions of 
the DMFT, such as the dynamical cluster approxima- 
tion (DCA) [14, 15]. Within these cluster approaches, 
characteristic difficulties have been encountered: The 



latest impressive work using the DCA by Maier et. al. 
performed a systematic cluster-size study of dSC in 
the 2D Hubbard model [14]. In clusters large enough 
(up to 26 sites) converged results point to a finite-T 
instability to dSC. Because of the QMC minus-sign 
problem, however, results were limited to U = At, 
where the typical energy separation in U and the mag- 
netic energy scale of the HTSC is not yet achieved. 
On the other hand, the VCA studies reviewed here, 
are clearly not yet converged with respect to the clus- 
ter size, as one can read off from Fig. 1 in Ref [6]. 
Cluster convergence is, at least in principle, also pos- 
sible within the VCA. With increasing cluster size 
longer-ranged correlations are included exactly. This, 
however, necessarily implies to use stochastic (QMC) 
methods as solvers for the cluster reference system. 
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